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Abstract 

The Shapiro Phase Shift on Pulsars (SPSP) influence at low , but 
detectable level the pulsars' (PSR) periods and their derivatives. This 
effect may be already written in PSRs' time derivative deviation and 
in particular it may explain the presence of a few negative pulse deriva- 
tives among 558 known PSRs. We describe in detail the phenomenon 
and predict that the SPSP may play a role in P deviations for just a 
few PSRs a year (~ 0.5 bevyr- 1 ) at aP ~ (l(T 14 1CT 15 ) ss" 1 
level , leading to a plausible explanation to the puzzling negative pulse 
derivative of few PSRs' periods. 

1 introduction 

The search for gravitational microlensing events of stellar sources produced 
by MAssive Compact Halo Objects (MACHOs) provides an interesting but 
approximated tool for the study of dark matter in our galaxy. Although sev- 
eral events of achromatic magnification have been already observed , they 
have not let the MACHOs' parameters (mass , transversal velocity ...) be 
known independently. Comparisons of the observed optical depths with 
those provided by theoretical galactic models have just suggested that there 
must be more matter toward the galactic centre than in a galactic halo , 
making more probable that our galaxy is a barred one. Therefore the real 
nature of the deflector that produces an observed event remains still unclear. 
Unfortunately the possibility of studying interference between the two beams 
produced by a deflector fails , due to the inchoerence of the stellar source. 
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Nevertheless the temporal phase shift produced by the deflector on " monochro- 
matic and coherent " PSRs' signals , the Shapiro effect , might provide an- 
other kind of investigation on the MACHO 's parameters. The extremely 
high precision the PSRs' period can be known , and the explicit dependence 
of the time delay on the deflector's mass , may offer a deeper tool for the 
measurement of both the deflector's mass and impact parameter. An event 
of temporal phase gravitational shift , due to the massive object crossing 
the ray's path between the PSR and the observer , produces the same phe- 
nomenon of displacing away and back the PSR on distances of the order 
of the Schwarzschild radius of the deflector. No Doppler shift (due to the 
gravitational pulling of a nearby crossing star) may lead to a similar forth- 
and-back phase shift , but it may only induce a one-time phase stretching. 
Doppler effects due to binary systems introduce a sinousoidal perturbation 
on residual PSRs' signal of the same order of magnitude but with a different 
time signature and a lower probability. 

A program of detection would require a continuous monitoring of the PSRs' 
period because every event is a unique one and leaves no trace ; however it 
may hold for a period long enough to be observed. 

We analyze the properties of this phenomenon by using the simple laws of 
gravitational lensing as found in literature. While the detectable gravitatio- 
nal magnification lensing needs millions of star sources (in order to fill up 
optical depths of ~ f0~ 7 ) , in our framework the Shapiro phase sensitivity 
(due to the logarithmic behaviour on the impact parameter) is much larger. 
This leads to larger detectable impact parameters , and so , larger optical 
depths (at least up to the needed r ~ f0~ 2 for the phenomenon to be ob- 
servable with hundreds of PSR sources). 

A similar effect between the two images produced by a microlensing event 
was already introduced || and it seems to be already observed on the PSR 
B0525+2f Q : it is our belief that the associated optical depth be too low 
(see note at proof). 

2 TIME DELAY VARIABILITY 

Consider the gravitational lens j|] due to a MACHO passing by the line- 
of-sight of a PSR with transversal velocity v±_ at a distance from the 
observer. Although two images are produced by the lensing of the PSR , 
the secondary image is observable only when the MACHO enters within the 
microlensing tube in general defined by the threshold parameter u = 1 
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|S, [U], |f| ; this means that the probability the secondary image be luminous 
enough is too low to take care of it , even if we assume at least ~ 500 
PSRs. Therefore from now on we shall consider only the primary image 
being observable. 

The formula approximated (see the exact one in appendix A) that gives the 
time delay of the primary image's lensed ray relative to the undeflected ray 



, very similar to the Shapiro one |11] , depends on both a geometrical and a 
gravitational term : 
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where u is the adimensional impact parameter of the MACHO from the 
line-of-sight of the PSR in units of the Einstein radius : u = b / Re (see 
figure 1) , where 



R, 



2r s 



D d D 



4.84 • 10 13 



M 



AKpc 



0.5 



0.5 



cm. (2) 



rs is the Schwarzschild radius of the MACHO , , D s are respectively 
the observer-deflector and observer-source distances and Dd s the deflector- 
source distance. In the Einstein radius formula the substitutions Xd = 
Dd/D s and 1 — Xd = Dd s /D s have been made. D s is expressed in charac- 
teristic PSRs' distance units. 

The geometrical term decreases very rapidly with the increase of u and it 
does not play any role as long as u > 1. The gravitational terms (a first one 
constant and a second one variable) account for most of the gravitational 
effect along the ray's path. The constant in the gravitational term under- 
lines the positive behaviour of the time delay. The logarithmic divergence 
for uRe S> Dd , Dds , i- e. for u — > oo , is fictious ; the time delay vanishes 
as ~ lf~ ' as shown in appendix A. As the formula shows , the time delay 
depends explicitly on the mass of the MACHO because of the presence of 
rs ; so it gives more transparent and accurate informations on the lensing 
dark object than the measurement of any gravitational lensing magnifica- 
tion does (oc 1/u) at small u. Equation (1) is valid nearly always (as long 
as u <C ^D s /(2rs) < 10 8 , for a characteristic solar mass) , namely for any 
realistic u < 10 3 . 

As the MACHO approaches to the line-of-sight of the PSR , the impact 
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parameter varies (in a first approximation) according to the law : 



u(t) = u rt 



1 + 
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where u m i n is the minimum impact parameter of the MACHO (figure 1) and 
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is a characteristic time. The largest value of t c is provided by a mass passing 
half-way between the observer and the PSR ( Xd = 0.5 ). A rough estimate of 
the integral phase residual time for our Shapiro shift in a 10 t c characteristic 
scale (an arbitrary cut-off time chosen for experimental data convenience) 
is : 

At (u min ) - At (10 

)-^(tT~) *"> . ( 5 ) 



M 



for large values of u m i n . Detailed residual time evolutions from equation (1) 
are described in figures (2) and (3) for deflecting masses M Q and 1O" 2 M . 
As we are allowing (on the contrary of the optical gravitational lensing 
of stars) large values of the impact parameter u we can neglect the first 
vanishing geometrical term in equation (1) and consider only the variable 
term in At grav . 

The speed of variation of the period P of the PSR is roughly the ratio of 
rs/c and the characteristic Shapiro phase shift event time , namely its time 
evolution is : 
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The considered i-b- varue f° r the present case , when the MACHO is half-way 
, has its minimum value. Therefore it defines a characteristic lower bound 
for residual PSRs' period variation. The function F(t) has a maximum and 
a symmetric minimum for the values =F tmax > where : 



1 + 



u„ 



(9) 



(the approximation being valid for large values of u m i n ) whose value (see 
equation 4) is nearly 4.6 y^jf^ yr • In figure (4) the function F(t) is 



yr 

y io< 

plotted for the values u m i n = 1 and u m i n > 10 . It can be seen that as the 
MACHO approaches to the line-of-sight of the PSR (negative values of t) , 
the PSR's period increases at a large rate , but , when u = u m i n the growth 
decreases to zero ; the process is time-symmetric. 

Using a double exponential disk model Up, the gravitational lensing optical 
depth can be evaluated. Choosing the galactic centre direction (/ = 0° , 6 = 0°) 
the calculation leads to the value 



2.33 • 10 u m i n 
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(being the integrated contribution of matter of the galactic bar component 
very low with respect to the galactic disk one between 4 and 8 Kpc from the 
galactic centre , it has been neglected). If this value is chosen as a reference 
one , both larger threshold parameters (of at least u-min = 

10 2 ) and the 

number of available PSRs ( > 500 ) can make the total optical depth be of 
the order of unity , i. e. a number 
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of gravitational host MACHOs can always be found at the same time. 
A rough estimate of the rate of "complete" Shapiro phase shift events is : 
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where in the last member the dependence on Xd has been eliminated by an 
exact volume integration. This rate estimate is based on known luminous 
lensing object galactic densities ; in presence of dark matter it may be one 
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or two orders of magnitude underestimated. 

A competitive gravitational delay or enhancement in residual PSRs' peri- 
ods may occur when an unbounded object (planet , star , black hole , ...) 
induce a gravitational deflection to the PSR (or the solar system) , leading 
to Doppler shifts (see appendix Al). Another equivalent shift occurs also 
for the presence of binary bounded objects near the PSR. Their influence 
induces a periodic sinousoidal residual signal in the PSR's period. Such an 
effect has been already observed at the level of intensity (/is over months) 
comparable to or below the one we are considering. However the signatures 
of such two different events are totally different (a step function increase or 
decrease , and a sinousoidal period residuals). The corresponding intensity 
and probability (related to the probability to find such an object in the sur- 
rounding volumes of the PSR or of the solar system) are extremely low with 
respect to the corresponding SPSP effect , whose probability is defined by 
much wider and longer cylindrical volume. 



3 CONCLUSIONS AND FIRST EXPERIMEN- 
TAL EVIDENCES 

The (SPSP) induces a unique characteristic secular variation on the PSRs' 
residual period : its signature leads to a pulse episodic bell-shaped time 
evolution of an order of tens of fis for one-solar-mass dark objects , well 
within the present detectability. It will never lead to an opposite sign event 
(first a decrease followed by an increase of the period). The pulse period 
variation could not be easily confused by a periodic Doppler shift , due to 
gravitational bounded planets , nor by rarest step increasing (or decreasing) 
Doppler shifts , due to an unbounded planet-star encounter near the PSR. 
The characteristic time residual At (equation 1) event , being dominated by 
a logarithmic function of the adimensional impact parameter u , has a wider 
cross section than gravitational lensing magnification and imply a larger 
probability to observe such an event even if at slow time scales. The charac- 
teristic scattering time is t c ~ 4.6 ( U JW ) (Jrj) \/ TT~ V r an d the cor- 



responding maximum amplitude in the period's derivative evolution 



P 



\F{t)\ P < 6.9-10" 14 (^gr) 1 (if^) yjj^ ss~ l which is above or within 
the present instrumental sensitivity and it is even comparable to the present 
characteristic PSRs' period derivatives. The characteristic time evolution 
shown in figure (4) makes the effect observable and well defined. However 
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Table 1: PSRs that manifest negative P. 
PSR P (ms) P (ss^ 1 ) Environment 



B 0021-72C 


5.7 


-4.0 


10" 


-17 


c 


B 1744-24A 


11.5 


- 1.9 


10' 


-20 


B , C 


B 1813-26 


59.3 


-3.0- 


10" 


16 




B 2127+11D 


4.8 


- 1.1 


io- 


-17 


c 


B 2127+11A 


110.6 


-2.1 


10- 


-17 


c 



In the last right column the notations mean : C = globular cluster , B = 
binary PSR 



in the present catalogues of PSRs the period derivative P has already an 
average positive derivative near 10 -14 s s -1 . This is related to the intrinsic 
PSR spin-down , due mainly to its luminosity radio emission. Therefore the 
Shapiro positive period derivative in earlier stages ( P > ) might be easily 
masked by common loss of angular momentum of the PSR. On the contrary a 
negative Shapiro period derivative near its maximum amplitude (\F\ ~ 0.5) 
in the second phase of the gravitational deflection (positive times , see figure 
4) might overcome the positive PSRs' period derivative leading to an excep- 
tional P. For u m i n = 10 2 we must expect nearly 0.5 ev yr~ l for M ~ M & 



at an intensity 



P 



< 6.9 -10 14 s s 1 and even more events (~ 5evyr 1 



for M ~ 10 2 M ) but at a lower level 



P 



< 6.9 • 10" 15 ss" 1 , almost 



on the edge of the characteristic intrinsic positive P. Therefore lower-mass 
deflectors ( M < 10 -2 M & ) will not easily perturbate the median intrinsic 
positive derivative. After a first inspection into the up-to-date catalogue 
of 558 PSRs E3] we do note the "anomalous" negative period derivatives 
shown in table 1. We note that the absolute P- values are quite small , but 
assuming an average positive P ~ 10 -14 s s^ 1 we are well within our rough 
predictions (0.5 4- 5 ev yr~ l ). We also note that three PSRs (first , second 
and fourth ones) are found in globular clusters , where is more probable to 
find a star near the line-of-sight. However the globular cluster optical depth , 
even if additional to the external galactic one , is not large enough to change 
significantly the SPSP rate probability. Finally we note that 85 PSRs has 
unmeasured P , possibly reflecting either intrinsic PSR stochastic behaviour 
or low level Shapiro phase shift noise. From figure (4) the maximum vari- 
ability occurs near the maxima or the minima , where P < . Any PSR 
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exhibiting the peculiarity P ~ , P < (or in the fast transient linear 
phase when P < , P ~ or , finally , in late stages when P ~ , JP > 0) 
must be carefully kept under control. Their i ? (t)-evolution may lead to the 
expected double or unique bounce from its maximum or from its minimum. 
A careful analysis and monitoring on P(t) evolution and deviation for all 
known PSRs (in particular the 85 whose P variability is unknown and the 
above-mentioned five negative ones) may offer a unique tool to observe dark 
matter MACHOs. Because of the large impact parameter , the SPSP may 
observe also diluited dark molecular clouds which might be spread over a 
wide radius ( 10 17 cm > R > 10 13 cm ) and could not be easily observed by 
gravitational lensing magnification. In conclusion we do believe that SPSP 
is already acting on PSRs' periods (see also note at proof) and it may better 
probe the mysterious nature of dark matter in both the galactic disk and 
halo. First evidences from a few low (10~ 16 -f- 10~ 17 ) s s~ l negative period 
derivative seem to favour light MACHO candidates of M < 1CT 1 M & with 
respect to large ones M > M@. Being the only one out of any cluster and 
any possible environmental disturbance , the PSR B1813-16 seems to be the 
best candidate suffering the SPSP process. 

A possible mimic of the SPSP may derive from the exchange of angular 
momentum from an accreting disk near a PSR , leading to 

6 

p =- 2 - 33 ' i(rl2 (£) MM' (i3) 

for characteristic PSRs with accreting mass fluxes of ~ 10 17 gs~ l . However 
such a flux would be characteristic of a pulsating X-ray source ; moreover 
the presence of such a disk (even of lower flux) near a non-binary object is 
quite unexpectable : neutron stars are assumed to evacuate external volumes 
after their super-nova birth. The alternative (and extremely rare) capture 
, destruction (by tidal forces) of a planet would introduce , by keplerian 
eccentricity , a recognizable periodic Doppler modulation. Hypothetical 
but unknown secular angular momentum exchange (between external and 
internal fluids of the PSR) might also lead to such a period deviations. 
The "artificial" fine-tuned accretion of a circular disk near the PSR remains , 
however , the unique alternative possibility which we were able to argue. 
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The gravitational time delay term is obtained by integrating along the ray's 
path the potential part of the effective refraction index introduced by the 
MACHO : At grav = - Jr J Udl = T -f J f 0. The calculation leaded to 
the subsequent more general formula : 



A GRAVITATIONAL TIME DELAY 
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being the coordinate of the primary image (in the 

uRe 



image plane having the origin towards the PSR , see figure 1) , b - 
the impact parameter and , therefore , (b + x+) = ^f- u + Vu 2 + 4 
As long as u — > oo , a\ ~ , 02 ~ J 

— ^— Therefore 
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as -u tends to infinity At 



grav 



4f and 

uD d 

fades to 



zero like ^ ^ (as it is expected). 

Equation (14) is surely too sophisticated for our purposes. A simpler for- 
mula can be deduced from equation (14). In our limit b, Re "C Dd,Dd s 

; therefore in this approximation only the third 



U <.r 



and fourth terms in equation (14) are important and , in both 03 and 04 
, the second term in the square parenthesis is neglectable with respect to 
unity. The resulting formula is equation (1). 



A.l THE DOPPLER SHIFT BY ENCOUNTERS OF NEARBY 
STARS 

The possible momentum exchange due to an unbounded nearby object (star 
, planet , black hole , ...) in gravitational deflection leads to a tiny Doppler 
shift on the PSR's period. Let us assume a scattering parameter b (not to be 
confused with the above-considered SPSP line-of-sight impact parameter) 
by a one-solar-mass object of Schwarzschild radius rs . The consequent 
orthogonal momentum transfer leads to an adimensional velocity exchange 
on the primordial PSR's motion is ||, |j : 

where (3@ is the adimensional relative velocity between the PSR and the ob- 
ject. The corresponding average period derivative (P) during the whole 

\ / Dp 

time AT of the gravitational deflection is 

(19) 
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This value which already assume a high and unprobable star density near a 
PSR ( > 10 5 pc~ 3 ) is negligible for most PSR's periods with respect to the 
characteristic SPSP delay. 

In addition to the orthogonal momentum transfer , the parallel momentum 
transfer , in analogy to the electromagnetic bremsstrahllung , can reproduce 
also a bell-like SPSP time delay. If we observe the PSR with the line of 
sight in the same direction of propagation of the nearby crossing object 
(MACHO , Juppiters ...) , because of the exchange of parallel momentum , 
the PSR is sent away and back by , roughly , a distance Al max m 9 2 M 2 in a 

characteristic time , producing a bell-like time delay whose maximum 

value is 

r s 1 

C P& 

Therefore the final figure will be a generic combination of the two effects 
(bell- like for parallel and step- like for orthogonal momentum transfer). We 
note that such bell-like PSR's period evolution may be sign-reversed because 
of the time-reversibility of the PSR-object scattering. In conclusion , out 
of a very rare geometrical configuration (in line of sight of the scattering 
direction) , the process should be in general distinguishable from our unique 
positive bell-like SPSP evolution. 



NOTE ADDED IN PROOF 

We acknowledge the referee for letting us notice a previous article [f?j where 
the SPSP have been applied to the residual period evolution of the PSR 
B0525+21 (D s = 2.3 Kpc , I » 183° , b = 6.9°). The authors found out 
a bell-shaped evolution which is similar to the SPSP effect. However , 
even if their formulation of the integral phase residual time agrees with our 



At{u) - At(u min ) « -t^ log 



(for large u) , they did not dis- 



cuss the optical depth for such a gravitational lensing event. From the 
fitting parameters of the deflecting mass ( M ~ 330 Mq ) and characteris- 
tic evolution time (t c ~ 0.32 yr) , equation (4) allows a minimum impact 
parameter of 

0.51 / f3± 



V lo- 



(21) 



As the most of the relative contribute £ to the optical depth is pro- 
vided when the deflector is nearly half-way the distance D s , we conclude 
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that u m in < 1 . The theoretical optical depth for the PSR B0520+21 
(nearly in the galactic anticentre direction) is r < 2.11 • 10~ 8 . Due to this 
extremely low probability as well as the absence of the unavoidable twin lens- 
ing image (u m i n < 1) , we were forced to reconsider a slightly more probable 
explanation of the Doppler shift encounter (see appendix Al). Contrary to 
our earliest considerations for the orthogonal momentum transfer (see equa- 
tions 18 , 19) , the parallel momentum transfer (equation 20) can reproduce 
a bell-like SPSP delay (with some additional step- like pattern). 
We have imagined a free-free gravitational scattering and we noted that , 
if the largest delay must be At» max 32 ms and if it is reached in nearly 



5yr , it might be produced by a mass of M « 10 3 (tfBw) M q with 



of the probability for such a nearby encounter should be « 5 • 10 -7 10~ 6 , 
respectively for realistic luminous or dark matter assumption ; still a small 
and puzzling value. 





However a rough estimate 
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x , image plane 




Figure 1: Relative position of source , MACHO , observer and primary 
image. As u — > oo , x + — > and the image corresponds to the source. v± 
accounts for the relative PSR-MACHO-observer orthogonal motion velocity. 
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